Fluctuations of the number of particles within a given volume in cold quantum gases 
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In ultracold gases many experiments use atom imaging as a basic observable. The resulting image 
is averaged over a number of realizations and mostly only this average is used. Only recently the 
noise has been measured to extract physical information. In the present paper we investigate the 
quantum noise arising in these gases at zero temperature. We restrict ourselves to the homogeneous 
situation and study the fluctuations in particle number found within a given volume in the gas, 
and more specifically inside a sphere of radius R. We show that zero-temperature fluctuations are 
not extensive and the leading term scales with sphere radius R as R^ In R (or In R) in three- (or 
one-) dimensional systems. We calculate systematically the next term beyond this leading order. 
We consider first the generic case of a compressible superfluid. Then we investigate the whole Bose- 
Einstein-condensation (BEC)-BCS crossover crossover, and in particular the limiting cases of the 
weakly interacting Bose gas and of the free Fermi gas. 

PACS numbers: 03.75.Hh, 42.50.Lc, 05.30.Fk 

I. INTRODUCTION 

In recent years major progress has been achieved in realizations of highly controlled low temperature quantum gases. 
One of the basic techniques, which permits the investigation of the properties of cold atomic clouds, is imaging of the 
condensate on a charge-coupled device camera. Usually the noise, present in each image of a particular realization 
of the cloud, is ignored and the overall density is fitted by a smooth function. This enables the extraction of the 
density profile and size of the cloud. Only recently a new generation of experiments has appeared [l| in which the 
noise itself is analyzed. It was shown f]\ that the study of the correlation functions in real space reveals correlations 
due to the presence of fermionic dimers in the BEC-BCS crossover problem. Naturally, proper control of these kinds 
of experiments requires a good theoretical understanding of the noise in these systems, and more specifically of the 
influence of quantum properties on the noise since these properties are the major interest of these systems. In the 
present paper we point out that these correlations at low temperatures indeed show interesting quantum properties 
and permit the extraction of additional information on the interactions in the system, for example, at unitarity. Since 
the systems of interest display the BEC-BCS crossover, we investigate in detail density fluctuations throughout this 
crossover. 

The problem of fluctuations was addressed in a number of theoretical studies 0, d, 0i [E| ■ Recently Belzig et al. 
studied fluctuations (^SN'^) = {{N — N)"^) in the BEC-BCS crossover at zero temperature using BCS theory. In 
their approach number fluctuations in a large enough volume grow linearly with the volume. It was predicted that, 
deep in the BCS regime, fluctuations are proportional to the value of the gap and are vanishingly small. However, 
because of their use of the simple BCS wave function, they did not take into account the collective mode, which 
k> as we will see plays a crucial role. Number fluctuations of a Bose gas in a sphere were studied with logarithmic 
accuracy in f3l|. Classical thermodynamic fluctuations are known to be extensive, i.e. additive with respect to the 
volume. On the contrary, at zero temperature fluctuations grow less rapidly, which is physically reasonable since the 
thermal source of fluctuations is no longer present. On the other hand there are still quantum fluctuations. It was 
found that these fluctuations at zero temperature are not additive and grow less rapidly than the volume. In the 
present paper we extend this work in two directions. First we go to the next significant order in the dependence of 
fluctuations on volume by performing the calculation with nonlogarithmic accuracy. Next we extend our study to the 
BEC-BCS crossover, which allows us in particular to cover at the same time the case of free bosons and free fermions. 
Fluctuations in systems with reduced dimensionality are analyzed by studying properties of a one-dimensional free 
Fermi gas. 

The paper is organized as follows. In Sec. |TT] we recall as an introduction to our calculations the basics of the 
fluctuation calculations. Then, in Sec. IIIIl we address briefly the case of a general compressible superfluid and recall 
that the dominant term is directly linked to long wavelength excitations, i.e. phonons. After this, in Sec. IIVI we 
proceed to a more precise calculation in the case of weakly interacting bosons and obtain explicitly the next term 
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beyond the dominant logarithmic one. Section |V] is devoted to the extension of this calculation to the BEC-BCS 
crossover, making use of the dynamical BCS model. Finally the case of the ideal Fermi gas is handled exactly in 
section IVTl both for a three-dimensional and a one-dimensional space. In the latter case we also comment on non zero 
temperature properties. 

Notice, that we will consider Bose systems only in the presence of interparticle interactions. Fluctuations in an 
ideal BEC gas are quite specific, due to the infinite compressibility of the system. A reader can find a comprehensive 
review of this subject in 



II. FLUCTUATIONS IN A SPHERE 

In a homogeneous system the average number of particles A'^ in a given volume within the system is merely fixed by 
the density and the size of the volume. In contrast, the particle number fluctuations depend on interaction properties 
of the gas. The strength of number fluctuations is quantitatively described by the dispersion (^SN'^) — {{N — N)^). 
These fluctuations are closely related to the density-density static correlation function ni^dri— r2|) = {6n{ri)Sn{r2)): 

{6N^)^n J j v{\v^-V2\) dv^dY2, (I) 

V V 

where n = {n{r)) denotes the average density in a homogeneous system, 5n{Y) = n(r) — n and V is the volume in 
which the fluctuations are measured. 

In homogeneous systems (to which we will restrict ourselves in the following) isotropy makes it natural to investigate 
number fluctuations in the situation where the volume is a sphere of radius R. Then the mean number of particles 
is A^ = |7ri?^n. The double integration in Eq. ([1]) can be partially performed resulting in a one-dimensional integral: 

2R 

(SN'^) = 47rn J i^{r)T{r)r^ dr, (2) 



where we have introduced the overlapping volume T(r) = Jy Jy dridr2(5(ri — r2 — r) of two spheres of radius R, with 
centers separated by r, which is given explicitly by 

r{r)^^-nrR^ + —. (3) 

The density-density static correlation function (pair distribution function) is related to the probability of finding two 
particles separated by a distance r. The function v(r) is defined in coordinate space. Its counterpart in momentum 
space is the static structure factor S{k), which we define precisely as 

Below we will use directly our knowledge of S{k) to calculate (JA^^). It is useful to introduce also the Fourier 
transform of T(r) which is immediately obtained from its above expression as T{k) = Jy Jy dridr2e'^^'^^~'^^^ and is 
given by 

T(fc) = [sin(fci?) - fci^cos(fci?)]^ (5) 

Then Eq. is replaced by 

n 



(SN^) = TT^ / S{k)T{k)e dk. (6) 







If the /-sum rule is exhausted by the collective mode branch of the excitations, the static structure factor can be 
found from the Feynnian relation 

hk'^ 

^ ^ 2mws(fc)' ^ ^ 

where m is the mass of an atom and ujs{k) is the excitation energy. 

The knowledge of the single branch excitation spectrum in a weakly interacting Bose gas allows us to calculate, 
using Eqs. ((T|)-(I7]), the dependence of the density fluctuations ((5A^^) on the strength of interactions and the radius of 
the sphere. 
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III. COMPRESSIBLE SUPERFLUID 



The long wavelength behavior of the static structure factor in any compressible superfluid is universal and cor- 
responds to the presence of sound waves with velocity c and linear excitation spectrum ujp}i{k) — kc. By using 
the Feynman relation, Eq. ([7]), which in the general case is valid in the k limit, we immediately find a linear 
dependence of the static structure factor S{k) = hk/{2mc). Substitution of S{k) into Eq. (|4]) gives the long range 
decay of the density-density correlation function v{r) oc l/r'* (see Sec. 87 in There are three terms in Eq. ([3]). 
The contribution of the constant term can be calculated using the identity Att v{r)r'^dr = S{k) \k=o = 0. Thus 

this contribution scales as J^^ i^{r)r^dr cx oc E?. The linear term causes a logarithmic ultraviolet 

divergence for r ^ in the integral of Eq. This divergence is due to a failure of the phononic description at 
large k and is avoided [H by truncation of the integral at some distance ^ of the order of the healing length. This 
results in a logarithmic dependence of (5N'^) on R, since 

i?2 J^-^^ j.(r)r3 dr cx i?^ j^'^ ^ = i?2 ln(i?/^). The last term 

is already converging and scales as R^ since R v{r)r^ dr cx r dr on R^. Hence the leading contribution to the 
integral comes from the term in T(r) linear in r. The final expression reads as Q 

(5iV2)^^^2i /^fV (8) 

with the coefficient C inside the logarithm undetermined, as it depends on the unspecified short-range behavior. 

The result Eq. ([5]) does not depend on statistics and is valid both in fermionic and bosonic superfluids. Thus at 
T = the increase of {SN"^) when R ^ oo follows the law R^ \nR. In contrast, at finite temperature according to the 
classical thermodynamic equation, {^SN^') = —kBTn{dV/dP)T oc i?"^, which is faster compared to Eq. ([5]). This was 
first noted in 

Notice that at finite, but low temperature Eq. ([S]) is valid at the distances R ^ fic/{kBT). For distances satisfying 
the opposite inequality, the density-density correlation function decreases exponentially (see [7], Problem 1 to Sec. 87) 
and the classical result for number fluctuations is valid. 

It is important to note that we consider the fluctuations of the total number of atoms in a given volume. Fluctuations 
in the number of atoms in the condensate and in the number of atoms out of condensate separately increase faster 
with R: (SNq) cx i?^ at T = and cx i?"' at finite T (for detailed discussions see [1| and [1]). 



IV. SYSTEM OF WEAKLY INTERACTING BOSONS 



One can easily calculate the coefficient C in a system of weakly interacting bosons with mass ms, where the average 
boson density ub satisfies nsa^ <C 1 (a being s-wave scattering length). These conditions correspond to the present 
experimental situation for trapped BEG gases, which are perfectly described by the Bogoliubov theory. The excitation 

spectrum in this theory has the form ujsik) = ^ (kc)^ + (5^^)^" and sound velocity is c = ^/AttWbo^ h/rriBa- The 
static structure factor is then found from the Feynman relation Eq. ([7]) which is exact in the framework of the 
Bogoliubov approximation 

Sb (k) = ^ ^ . (9) 

The small-momentum part fc ^ describes phonons 6*3 (fc) ~ hk/{2mBc), while the large-momentum part fc — > cx3 
corresponds to free particles with energy given by = h?k^ /2mB and static structure factor ^'^(fc) 1. 

By substituting Eq. @ into Eqs. ©-(HI) we obtain the density fluctuations in a weakly interacting Bose gas (details 
of the calculation are given in Appendix |^ : 

(SN^-)^^R^l^(8e^-^^]^ (10) 
niBC \ " / 

where 7 = 0.577... is Euler's constant. We see explicitly that knowledge of the short wavelength physics is important 
for the calculation of the coefficient inside the logarithm. In our case this is done by assuming the Bogoliubov 
excitation spectrum. 

Notice that our theory is valid for a uniform gas. However, it can be applied to a small volume inside a trapped 
BEG gas in the Thomas- Fermi regime under the condition i? <C i?o, where Rq is the condensate radius. 
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V. SUPERFLUID SYSTEM OF TWO-COMPONENT FERMIONS 



We now consider the calculation for a system with two fermionic species in the BEC-BCS crossover within the 
dynamical BCS model. As we have mentioned it is convenient to use the numerical results which have already been 
obtained for S{k) in the study of molecular signatures in this crossover Q. Since we now deal with two kinds of 
fermions with the same mass mp and equal populations, say spin-up and spin-down fermions, we first have to take 
this into account in the general formulation of our problem outlined in Sec. |TT1 Following the same steps it is easy to 
see that the total number fluctuation is still essentially given by Eq. ([6]): 



oo 

(SN^) = ^ y S{k)T{k)e dk 



(11) 



provided we use the total density fitot = 2n, where we keep the notation fi = ^ for the single species average 
density. Here S{k) is the static structure factor defined by S{k) = 'S'ii(fc) + 5|j(fc), which has been calculated in 
Ref. [9] and is again related to the sound velocity c in the limit /c ^ by S{k) = hk/ {2mpc). Naturally we now have 
SN = SN-^ +SNi. 

Although S{k) is known numerically, brute force numerical calculation is not a good choice and it is better to 
extract analytically the R^lnR behavior which comes from the low k domain, as we have seen quite generally in 
Sec, mil This is most conveniently achieved by introducing a simple model So{k) for the structure factor which can 
be handled analytically. Then the difference with the actual S{k) is handled numerically. We take for So{k) a model 
which behaves as iS'(fc) in the two limits /c — > and k — > oo, in order that the difference S{k) — So(k) goes to zero in 
these two limits. This simplifies both the general analysis and the numerical calculation. We choose: 



5o(fc) = ak, 
So{k) - 1, 



fc<i 

Q 

fc> i 



(12) 



where we take the coefficient a = h/{2mpc) in order to recover the proper behavior for k 
analytical calculation of ((57V^) for this model are given in Appendix lAl The result is 



0. The details of the 



(SN'^) ^ = SnaR^ In 

\ / mod 



R 

In - + 1.270. 

a 



(13) 



Then we are left with the calculation of 

' l/a 



dk ^^^\. [sin(fci?) - kRcos{kR)f 



fc4 



^j^S{k}l jgjj^^^^^ _ kJicoa{kR)f 



l/a 



.(14) 



In both terms the dominant contribution in the limit i? — > oo comes from the kRcos{kR) term, and moreover, writing 
cos'^{kR) = (1/2)(1 -t- cos(2/ci?)), we see that only the 1/2 is left at dominant order. In other words we may replace 
cos^ (kR) by its average 1/2. This leads to 



mod 



l/c 



dk 



S{k) — ak 
A?2 



dk 



S{k) - 1 
fc2 



l/a 



(15) 



We see in particular that the first integration no longer gives rise to a divergence in the limit fc — > 0, which is present 
when one has S{k) instead of S{k) — ak in the numerator and is responsible for the R^ Ini? contribution. Hence the 
integrations are regular in both the limits fc — > and fc — > oo and are easily performed numerically once the numerical 
values of S{k) are known. The result is presented in Fig. [TJ We have written the number fluctuation in the general 
form (SN"^) = 8naR^ [\ii{R/a) + C], where a = h/{2mFc) and we have displayed the coefficient C as a function of 
l/{kpa). We have also plotted in Fig. [21 from this asymptotic result, the behavior of the number fluctuation as a 
function of the radius of the sphere i?, for various values of l/{kpa). Our extension of this figure toward fairly low 
values of kpR is justified by our finding in the case of the free Fermi gas, considered in the next section. 

It is naturally of interest to consider the BEG limit \/{kpa) -\-oo and to compare the numerical result with the 
analytical one found in Sec. |TVl First the number SN considered in Sec. [IV] corresponds to fluctuations of bosons. 
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FIG. 1: Two-component Fermi gas in the BCS-BEC crossover, where the coefficient C in the logarithm is defined as (SN^) — 
8naR^[ln{R/a) + C] with a = h/{2mFc). 




FIG. 2: (Color online) Density fiuctuations in two-component Fermi gas in the BCS-BEC crossover: dashed line, l/kpa = —2 
(BCS regime); solid line, l/kpa = (unitary regime); and dash-dotted line, l/kpa — 1 (BEC regime). 

Since in the BEC limit the bosons we end up with contain two fermions, we have to multiply the result of Sec. IIVI bv 
four in order to obtain a fermion fluctuation which can be compared with Eq. (jlOp . Hence we expect to find in the 
BEC hmit 

/SN'X^^^R^l^fiee'^-^^) (16) 
mpc \ / 

since with our convention we have hs = n, while tub — 2mp. Taking into account a — h/(2mpc) we see that the 
cocfiicicnts in front of Ini? agree as they should. Next from Eq. (|16p we should find C = 31n2-|-7— 1 = 1.657... 
for the coefficient C defined above, while we see from Fig. [2]that for 1/kpa = 4 the result is C = 1.625.... Since we 
are not yet in the strong BEC limit 1/kpa oo, the agreement is excellent. It is worth noting that a very good 
agreement with the BEC limit is already found starting from 1/kpa ~ 1. 

It would also seem natural to investigate in the same spirit the BCS limit 1/kpa — > — oo. However, it is well- 
known that this limit is somewhat singular. Indeed, when 1/kpa goes toward — oo, S{k) becomes identical to the 
structure factor of the free Fermi gas; but in the small k range this occurs in a singular way. For a free Fermi gas, 
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FIG. 3: Qualitative behavior of S{k) at low k in the limit of large and negative l/kpa. In the limit fc ^ 0, S{k) = ^/3k/2kF 
(dotted line). For k ~ 1/^, where ^ is the coherence length, this behavior switches to the free Fermi gas result S{k) = Sk/Akp 
(dashed line). 



S{k) ~ 3k/Akp for k ^ 0, that is a = 3/Akp, but the sound velocity we will find for large negative 1/kpa is given 
by c = \/ {n/mp)dfi/dn = Ukp/mF^/i which means that the relation a = hjlmpc does not apply. Hence, we can 
not recover the free Fermi gas in the limit l/kpa — >■ — oo since we will already have a disagreement at the level of the 
coefficient of the R^lnR term. The basic physical reason is that the above expression for the (first) sound velocity 
implies naturally thermodynamic equilibrium, while in the non-interacting gas the relaxation time is infinite which 
prevents the existence of first sound. On the other hand it is known that in the superfluid state this singular behavior 
disappears and first sound exists with a velocity related to compressibility and given by the above general formula. 
For large (but finite) negative l/kpa, S{k) has the following anomalous behavior, which is displayed qualitatively in 
Fig- [3 For A: — > one has indeed S{k) ~ fik/2mpc — \/3fc/2fci?. However this behavior holds only in the small range 
Q ^ 1/C ~ AmF/hkF, where ^ is the coherence length and A the BCS gap. Beyond this range one recovers the free 
Fermi behavior S{k) ~ Sk/Akp; but, since A — > when l/kpa — oo, all this structure for S{k) occurs in a range 
of vanishing size in this limit. If we insist on extending our above calculations in this range, the change of slope of 
S{k) at small k we have just described will produce a divergent result for our coefhcient C in the limit l/kpa — oo. 
Another way to state this result is to say that the two limits 1/kpa — > — oo and i? — > oo do not commute. One finds 
the free Fermi gas limit if one keeps R finite and first let l/kpa — oo, but not when these limits are exchanged, as 
is clear by considering the coefficient of the R^ In R term. 

After having made this point we calculate (JiV^) nevertheless for completeness in the following section for the free 
Fermi gas. 



VI. IDEAL FERMI GAS 



In the case of ideal Fermi gas it is possible to perform a fully explicit calculation of SN'^. We can then compare the 
logarithmic-accuracy expression with this exact result. In an ideal Fermi gas there are no correlations between atoms of 

different spin orientations, and hence nv^i{r) = ((rif (r) — n)(n|(0) — n)) = and (^SN^^'j = ((iVf — iV|)(A^x — iV^)) = 

0. The same-spin density-density correlation function of an ideal Fermi gas at zero temperature is known 1(1 explicitly 
and manifests decaying oscillations [l^ with a period defined by the Fermi momentum kp — (67r^n)^/^: 

-n{r)=5\r)-^^[-^-coskpr) , (17) 

describing the Fermi hole due to Pauli exclusion, with the free particle behavior being recovered in the r —^ range. 
Then density fluctuations inside a sphere are given by Eq. ([2]). The integration can be performed exactly in terms 
of special functions and the full expression is given in Appendix [B] by formula (jB4p . Here we are interested in the 
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FIG. 4: (Color online) Number fluctuations in an ideal Fermi gas. Solid green line - exact evaluation of the integral Eq. HB4|) 
and dashed red line - leading logarithmic contribution Eq. (|18p 



long-range asymptotic. The leading term is 

(SN^A = -^^^-^ In lie^-ikpR] + OikpR). (18) 

The comparison of the asymptotic behavior Eq. (fTS]) with the exact result (|B4I) is shown in Fig. ID An important 
observation is that this long-range asymptotic works very well for distances as small as kpR w 1. The coefficient 
under the logarithm is 4e'''~2 ~ 4.32. 

Of course, the general logarithmic accuracy expression Eq. ([5]) coincides with the direct calculation Eq. (|18p once 
expressed in appropriate units. We note that for an ideal Fermi gas one should use c = 2cF/3m (where cp = hkp I'm 
is the Fermi velocity) as the expression for the speed of sound, as this quantity enters in the long-range asymptotics of 
the density-density correlation function and corresponds to the slope of the linear behavior for small momenta in the 
static structure factor. This velocity is different from the speed calculated through the compressibility m(? = nd^/dn, 
which leads instead to c = cplsf^. The difference is quite natural because there are no sound- like excitations in an 
ideal Fermi gas. This difference is absent in bosonic systems and the system of ideal low-dimensional fermions. 

At unitarity the only relevant length scale is again fixed by the Fermi momentum (or density), while the chemical 
potential is obtained from the non interacting one by multiplication by the coefficient ^ ~ 0.42(1) p^. Accordingly 
the speed of sound is given by c = \J\Ukp j \fZm. We note that the value of ^ can be extracted from experiment by a 
precise measurement of the total fluctuations, since at unitarity we expect 

(^ln[Cfc^i?]. (19) 

Analogous results can be obtained for an ideal one-dimensional Fermi gas in a box of size 2i?, as considered in 
detail in Appendix [Cl where the nonzero temperature case is also considered. The fluctuations are given by 

{SN'^) = z/(r)T(r) dr, (20) 



where T(r) = 2R — r and n is the particle density by unit length. For an ideal Fermi gas the density-density 
correlation function at zero temperature is i/(r) = S{r) — nsin'^{kFr)/{kFf')'^, where kp = 7rn is the one-dimensional 
Fermi momentum. The leading term is logarithmic and, with inclusion of the next order correction, the result is found 
to be: 

{SN^)[l^ = ^ HAe-'+^kpR). (21) 



The complete expression for (^SN^') is given by formula (jC6p . 
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The low temperature expression for the density-density correlation function i/(r) is obtained in Appendix [C] and 
is given by Eq. (|C8p . Number fluctuations at finite temperature are obtained from the integral ([^0)1 . Numerical 
integration shows that at finite temperature the additivity is restored and number fluctuations are proportional to R. 

Notice, in conclusion of this section, that the asymptotic law (JTV^) ~ R-^\nR is valid also for an arbitrary non- 
superfluid Fermi liquid, because one has always S{k) ~ fc at fc — > 0. The prefactor in this equation can be calculated 
in terms of the /-function of the liquid by solving the Landau kinetic equation (see Sec. 91 in Q)- However the 
coefficient inside the logarithm depends on the high momentum behavior and cannot be calculated in a general form. 



VII. CONCLUSION 



In conclusion we have investigated in this paper the zero-point fluctuations of the particle number in a macroscopi- 
cally large sphere of radius R for different types of superfluids. In all three-dimensional cases at zero temperature we 
obtain an anomalous dependence ((57V^) ~ R' \nR. The prefactor in this equation depends on the system, but can be 
calculated in general form, while the coefficient inside the logarithm depends on the high momentum behavior of the 
system under consideration. We have calculated this last coefficient explicitly throughout the BEC-BCS crossover, 
and in particular for the limiting case of the weakly interacting Bose gas. We have also considered the somewhat sin- 
gular limit of the non-interacting Fermi gas. Our obtained results can be relevant for noise measurement experiments 
in cold gases. 

We thank S. Stringari for useful discussions and Yu. M. Tzypeniuk for pointing out Ref. [4]. We thank B. Jackson 
for reading the manuscript. G.E.A. acknowledges support from MEG (Spain). 



APPENDIX A: DERIVATION FOR THE MODEL STRUCTURE FACTOR AND FOR THE WEAKLY 

INTERACTING BOSE GAS 

From Eqs. (O and we have to calculate, with the change of variable x = kR 

oo 

J [sin(fci?) - kRcosikR)f = aR^h + R^h, (Al) 





where 

R/a 

f 1 2 

Ji = / [sinx — xcosa;] , (A2) 

J x-^ 



CX3 

f 1 2 

I2 ~ / dx—[smx — xcoax] . (A3) 

R/a 

Since we are interested in the limit i? cx), it is clear that in I2 only the contribution from a; cos a; is left at the 
dominant order in R we are interested in. Moreover, we can, for the same reason, replace cos^ x by its average 1/2, 
which gives to dominant order 

(A4) 

It is in principle more complicated to handle Ii since the dominant In R term comes from the contribution of x cos x at 
the upper bound, the other terms becoming negligible in this limit, while for a; ^ all the terms need to be retained 
in order to avoid singularities. Nevertheless this can be conveniently done by integrating by parts: 

R/a 

1 2-^1 /" J sina;(sina; — cccosx) 



/i = — cos I dx 

2 a J 




2R/a 

9 + 9 / '^^ 7 = -ln— +7-C1 — -1. A5 

2 2 t 2 a a 
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Here we have omitted in the first step terms going to zero when R ^ oo since these terms are indeed neghgible at 
the order we are interested in. The last step introducing the cosine integral Ci(a:) can be found in Ref. JJJ. Since 
Ci(a;) — > when x — > oo we obtain by gathering all the terms 



aR^h+R'h^laR' In^. 

2 a 



This gives for the number fluctuation 



(SN^) = 8naR^ In = 8naR^ 

^ ' a 



In - + f.270 

a 



(A6) 



(A7) 



as indicated in the text. 

To handle the bosonic case where 



Ssik) = 



hk 



1 + ( nk ) 



(A8) 



we proceed by the difference with the above model, where we take naturally a — h/(2mBc), just as indicated in 
Sec. |Vl We are left with the calculation of: 



mod 



l/c 



dk 



S{k)-ak f „ S(k)-l 
' ' dk- 



fc2 



1/a 



With the change of variable k = 2niBCx/h = x/a we are left to evaluate the simple integrals: 



(A9) 



dx{- 



Vl + X 



-) =ln2-ln(l + \/2) 



dx{- 



^ - 4) = + \/2) - 1, 



Vl + a;2 x^ 



(AlO) 
(All) 



which leads to the result indicated in the text: 



Ae'^-^R 



a 



(A12) 



APPENDIX B: DERIVATION IN AN IDEAL FERMI GAS 

The same-spin density-density correlation function is given by Eq. ()17|) . from which the number fluctuations 

is obtained from Eqs. ([2]) and ([3|): 

2R 



AirnR^ 4 f dx f 3x \ , . „ . ^ „ ^^7 

' ^ / — ( "1 + — " Ye ) [sHkpRx) - {kFRx)cos{kFRx)Y 



(Bl) 



where the first term comes from the (r) contribution in (r) in Eq. (jl7p and we have just made the change r = Rx. 
Taking into account 



oo 

/ — [sin(fc_Fi?a:) - [kpRx) cos{kFRx)f = jikpRf, 

J X D 



(B2) 
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leading to a contribution which cancels the first term of Eq. (|B1[) . we find: 

2 



{^^^^) ^ ^ i — - —j [smikpRx) - ikFRx)cos{kFRx)f 



oo 

i f dx 

+ — - / —\sm{kFRx) — {kFRx)cos(kpRx)] . 



(B3) 



This expression can also be obtained directly from Eq. together with the expression of the structure factor for 
the ideal Fermi gas S{k) = (3fc/4fcj^) - k'^/{l&k%) for k < 2kF and S{k) = 1 for fc > 2kF- 
Setting R — AkpR, we obtain the following explicit expression [llj: 



(SN') = 



R^ 


'AiR) - \ 


i?2 


327r2 


^ 1447r2 




1 

"^2887r2 


[iRsinR- 



ttR - , 

— - RSi{R) - cosR 



(B4) 



where Si(a;) = J dt and Ci(a;) — ~ J dt are sine and cosine integral functions and A{R) = Ini? + 7 — Ci(i?), 

_ ^ _ 

with the dominant behavior In i? + 7 for large R. 

The leading term in a sphere of large radius R is derived from the first term and is indeed 



{^^11) = ^^^^^ 



Ae^-^kpR 



OikpR). 



(B5) 



APPENDIX C: DENSITY FLUCTUATIONS OF ONE-DIMENSIONAL IDEAL FERMI GAS AT A 

FINITE TEMPERATURE 

It is of interest to discuss the effect of temperature. The one-dimensional ideal Fermi gas provides a simple example 
where this can be done. A nice expression for i/(r) at low temperature can be derived, but still one has to calculate 
an integral involving i'{r) in order to obtain ^JiV^^. 

The density-density correlation function is written in second quantization as 



ni/(x,x') = (n(x)n(x')) - = (*^(x)*(x)*^(x')*(x')) - 



or, with ^'(x) — X^k^ "^k/Vv , we have in terms of plane waves creation and annihilation operators: 



ni^(x, x') 



1 

y2 



E 



ki,k2,k3,k4 



(CI) 



(C2) 



Using Wick's theorem we get (^al^ak^al^a]^^'^ = (^ki'^ks) (aLak^ 
7iki(5ki,k4(l~'^k2'^k2.k3)- In a homogeneous system the correlator (jC2l) depends only on the relative distance r = x— x': 



"ki°k4) (oksaka) = "ki^ka'^ki.ka'^ka.ki 



n;/(r) = n6^ {r) 



Uke 



Uke 



d^k 



(2^) 



D 



(C3) 



where D is the space dimension. For free fermions Uk = [e^*^* ''^/^ -|- l] .At low temperature the chemical potential 

can be well- approximated by the zero temperature value fi — . In a one-dimensional system one has to evaluate 
the following integral: 



cos kr 



e 2™t 



dk 
'2^' 



(C4) 
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Let us first consider the T = case. We find easily, as indicated in the text 

i.(r) = 5(r) - n— — (C5) 

The strength of fluctuation can then be obtained exactly by calculations analogous to the ones performed in Ap- 
pendix |B] We find, with the same notation R — 4kpR as in Appendix IB] 

TT^ (SN^) =lni? + 7 + l + ^- i? Si(i?) - Ci(i?) - cos(i?), (C6) 

where Ci(a;) and Si(x) are cosine and sine integral functions respectively. The first three terms in Eq. (|C6p provide the 
dominant contribution and are reported in Sec. IVII The contribution of remaining terms scales like — cos{R)/R'^ for 
R^ 1. Indeed, there are canceling contributions as can be seen from asymptotic expansion Si{R) = tt/2~ cos{R)/ R — 
sm{R)/R'^ + 2cos(i?)/i?3 + ^nd Ci{R) = sm{R)/R - cos{R)/R^ + 0{R-^). 

We turn now to the evaluation of v[r) at low temperature. As can be seen from integration by parts, the main 
contribution to the integral \CA\ comes from k ^ kp. We proceed with the integration in a dimcnsionlcss variable 
X = h^(k^ — kp)/2mT « h?'kp{k — kF)/mT and extend the lower limit of integration in Eq. (IC4p to minus infinity, 
thus adding exponentially small contributions. 

The resulting integral can be integrated by parts. 



oc 



^ _ fcA f cos{kFr + kxrx) _ 1 f sm{kpr + kxrx) _ kx sin kpr ^^^^ 
TT J + 1 "^r J (e^ + l)(e^^ + 1) sinhTrfcAr' 

— oo — oo 

where kx — is a characteristic thermal momentum of the problem. Now we can write the density-density 

correlator: 

I . ... _ / TxkxsmkpT y 
v{r) = 8{r) - n - — . . (C8) 

\kp smhTTK^r / 

This shows that the density-density correlation function in a one-dimensional Fermi gas at finite temperature has an 
exponential decay at large r (like in the 3D case at finite temperature), in contrast with the T = result, Eq. (|C5p . 
exhibiting a power law decay recovered from Eq. (|C8p in the limit of zero temperature kx 0. 
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